The role of Mie scattering in the seeding of matter-wave superradiance 
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Matter-wave superradiance is based on the interplay between ultracold atoms coherently organized 
in momentum space and a backscattered wave. Here, we show that this mechanism may be triggered 
by Mie scattering from the atomic cloud. We show how the laser light populates the modes of the 
cloud, and thus imprints a phase gradient on the excited atomic dipoles. The interference with the 
atoms in the ground state results in a grating that in turn generates coherent emission, contributing 
to the backward light wave onset. The atomic recoil 'halos' created by the Mie-scattered light exhibit 
a strong anisotropy, in contrast to single-atom scattering. 
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Matter wave superradiance (MWSR) [l[ and collective 
atomic recoil lasing (CARL) are light-induced insta- 
bilities of the density distribution in atomic clouds. More 
precisely, they are due to correlations between successive 
scattering events mediated by long-lived coherences in 
the motional state of an (ultracold) atomic cloud or in 
the light field of an optical resonator {ij. Despite con- 
siderable theoretical efforts having been devoted to the 
dynamics of MWSR [B-Q, open questions still remain. 
One of them concerns the seeding mechanism which is 
able to start the MWSR instability even in the presence 
of losses. Thermal and quantum fluctuations will natu- 
rally contribute. However, in this paper we point out the 
particular role of Mie scattering, which turns out to be 
important at the onset of MWSR. Indeed, Mie scattering 
is active before any instabilities have developed, and it 
induces a phase correlation between the atomic dipoles 
that favors the build up of an instability. 

The below-threshold dynamics and the seeding of 
matter-wave superradiance are interesting problems. As 
long as we consider the atomic cloud as a homogeneous 
entity, e.g., a Bose-Einstein condensate (BEC) in the 
mean field description, no scattering should occur at 
all. Theoretical models which describe BEC's as mat- 
ter waves without fluctuations thus fail to explain how 
MWSR is initiated in the absence of a seeding wave §. 
In contrast, recent work has shown 0, [l^] how atomic 
coarse-graining, density fluctuations and Mie scattering 
from finite-sized clouds can even influence the scattering 
of a single photon. Also optical cavities may strongly 
affect the scattering by shaping the angular distribution 
of the density of modes that are capable of receiving the 
scattered photons ll|. These processes have a decisive 



Here we show that Mie scattering, caused by the finite 
size of the atomic cloud, favors the formation of a matter 
wave-dipole grating. This has already been suspected in 
12| . Indeed, prior to any significant motion of the atoms, 
their dipoles collectively order, which in turn leads to 
coherent emission. 



impact on the mode competition preceding the exponen- 
tial instability, and thereby on the instability itself. 



In previous papers [1 31 - 41 61 ] , we discussed the impact of 
atomic coarse-graining and finite scattering volumes on 
the radiation pressure force which acts on the cloud's 
center of mass. Here, we investigate the momentum 
distribution following the cooperative scattering of the 
laser light by a BEC. Our theoretical model describes the 
atomic cloud as a macroscopic matter wave that is homo- 
geneously distributed within a sphere, i.e., the atoms are 
considered to be strongly delocalizcd and density fluctu- 
ations are neglected. We find that the momentum distri- 
bution of the atoms adopts the shape of a recoil halo, 
very similar to the ones observed experimentally in time- 
of-flight images of BEC's. The halo indicates the direc- 
tions into which the atoms are preferentially scattered 
before the density distribution is noticeably modified. In 
particular, it exhibits a pronounced peak at 2hk. This 
corresponds to an increased backscattering of light that 
acts as a seed for MWSR. 

We consider an homogeneous spherical cloud of two- 
level atoms illuminated by a laser of wavevector ko = 
koz, Rabi frequency fio = dE^/h (where d is the electric 
dipole matrix element and Eq the laser electric field) and 
detuned from the atomic transition by Ao. The atomic 
cloud is described as a bosonic ensemble of N two-level 
atoms with field operator ^(r,t) = $ g (r,t) + $ e (r,t) (g 
for the ground state, e for the excited one). We treat the 
condensate as an ideal gas and consider the scattering 
between matter waves and optical waves, but we neglect 
nonlinearities due to atom-atom interaction. In second 
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quantization, the interaction between the atoms and light 
is described by the Hamiltonian fl3, 17|: 



i/j e (r,t)/ip g o(r), with |/3(r)| 2 the probability for an atom 
to be excited, one obtains [28| : 



H(t) = — I dr %l(r,t)$ g (T,t)e- 



(1) 



iA fc t+ik-r 



h.c. 



where gk = d[u}k/(2HeoV v )] 1 ' 2 and V v is the quantiza- 
tion volume. In Eq.([T]), the first (second) line describes 
the absorption and emission of a pump mode fio (& vac- 
uum mode Sk)- Doppler effects are neglected. Replac- 
ing ^ e (r,t) — > ^ e ( r ) t)e lA °' induces an energy shift of 
— Ao|\l/ e | 2 in the Hamiltonian, from which the Heisen- 
bcrg equations may be derived: 
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i(A -A fc )t+ik-: 
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For large atom numbers and far detuning from the atomic 
transition frequency, one can neglect quantum fluctu- 
ations and treat the operators as c-numbers — > ij), 
flk — > ak)- Eq.(j4]) is integrated over time as 



-«A i (3) 
~ig^e- l{A "- Ak)t f dr* t „(r,t)* e (r,t)e- lk r . (4) 



dt / e -i(A -A k )f / dr^(r,t')Ve(r,t)e- 



ikr 



(5) 

and inserted into Eq.((3]). Furthermore, we switch to a 
continuous-mode description — > (V^/(2tt) 3 ) J dk and 
obtain: 

dt 



x y dk^e ik -( r - r ') y dt'e^-^^^'V^r'.^V-eCr'.t'). 

We consider time scales over which the atomic density 
does not significantly change, i.e., ip g (r,t) ~ ip g o(r) in 
Eq.©, with po(r) = |?/> g o( r )| 2 the initial density of the 
cloud. Using the Markov approximation, i.e., the pho- 
ton time of flight through the cloud is much shorter 
than the atomic decay time, the last integral in Eq.© 
is replaced by S(k — fco)V'go( r ')' ! / ; e(r', t)/c. With the as- 
sumption that all the electromagnetic modes are equally 
present in the system (g^ w gk ) and by keeping rotating- 
wave-approximation terms, one can show that (l8j 



dk? k e 



2 ik-d 
k r 



d £/ e i(A -A fc )(i-t') _ 



2ifc |d| 
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where T = V v g\ kg/nc = d 2 k^ / (2irheo) is the atomic 
decay rate. For the normalized excitation field /3(r, t) = 



gfet) A A r 



/3(r,t) 



,f2 



e ik -r 



, , exp(ifc |r - r' 
dr p Q (r )- 
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-/3(r',t), (8) 



which recovers the model of cooperative scattering of a 
plane- wave introduced in [l3j . 

In the steady-state regime, Eq.(j8]) can be used to de- 
scribe the scattering of light from a dielectric medium 
with an index of refraction 11611 



i-4^ r/fc 3 (2A + ir). 



(9) 



In this regime of linear optics, the field j3 is directly pro- 
portional to the electric field, and the ratio between both 
is the susceptibility. Hence, the excitation pattern inside 
the cloud can be calculated analogously to Mie's the- 
ory [3 H3- The polarization amplitude /3 inside the 
cloud is decomposed into modes, thereafter labeled n, 
which arc elementary solution of the Hclmholtz equation: 



(A + mj!*g)0 = 0. 



(10) 



In spherical homogeneous distributions (with coordinates 
r(r, 9, (p), where 9 is the angle with respect to the z axis) 
these solutions are of the form j n (m c kor)P n (cos 9) in- 
side the cloud, with j n the spherical Bcsscl function and 
P n the Legendre polynomials. Then the amplitude (3 n 
of each mode are calculated according to the method de- 
scribed in Ref. [lit ], giving for an uniform spherical cloud: 



/3(r) 



Y j (2n + l)i n Mn{m c k r)P n {cos6), (11) 



n=0 



zA ^ e (r, t) - i^ e *o-^ fl (r, f) - ^(r, t) / dr'(6) p n 



where the coefficients /3„ are given by 
j n (k R) 



(25 + i)j n (m c k R) + i\nh£\k R) 



(12) 



with R the radius of the sphere and hn the spherical 
Hankel function. In general, the radiation is accurately 
described by considering a number of modes on the order 
of the size of the system k^R [20| • Note though that such 
a coarse-grained approach to describe the cloud cannot 
capture features associated with disorder (e.g., Anderson 
localization). 

As can be deduced from Eq. (|10|) . the laser imprints 
a phase with wavevector m c ko to the excited state ip e , 
which interferes with the ground state tp g that has a con- 
stant phase. The resultant grating is at the origin of 
coherent emission and of the backscattering wave that 
acts as a seed for MWSR. Indeed, the field radiated by 
the atoms inside the cloud in a direction u(#, (f>) is pro- 
portional to the structure factor: 



Sc (u) = i J P {v)p{vy 



-im c feou-r 



dr, 



(13) 
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with the index of refraction m c given above. Hence, the 
coherent emission by the cloud is a direct consequence 
of the periodic excitation field and the resulting grating. 
Note that we do not perform an adiabatic elimination of 
the excited state [2l|; we rather extract the scattering 
pattern from this state. 

The scattered light is difficult to observe directly. 
However, the radiation pattern is also present in the 
momentum distribution of the atoms, which can be 
easily recorded by time-of-flight imaging. Different 
from an iV-body model, the quantum matter field ap- 
proach yields the momentum distribution simply through 
the Fourier transform of the matter field, ~ip(p) = 
(2ttH)~ 3 J ip{r)e~ %p ' r / h dr. Hence, for an homogeneous 
cloud, the momentum distribution of the excited state 
is directly proportional to the structure factor: V'e(p) °t 
s c (p/ft), with p = m c fik and k = kou. Using Eq. (fTTj) , it 
can be deduced that the momentum wavefunction rp e for 
an uniform sphere of radius R reads: 



V'e(p) 



r(2*n)» t 



^(27i + l)/3„7„(p)P„(cos0), (14) 



with 7 n (p) = 47r f* 1 r 2 j n (m c k r)j n (pr /ft) dr. For a very 
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FIG. 1. Momentum distribution of the excited state D e = 
\ijj e \ 2 in logarithmic scale. The emission pattern |s c (k)| 2 con- 
centrates around a circle with radius p — m c hko (black line). 
Simulations for a laser detuning Ao = — 3GHz and a cloud (a) 
of size koR = 20, atom number N = 1.15 x 10 6 and refractive 
index m c = 1.2, (b) of size koR — 20, atom number N = 100 
and refractive index m c = 1 + 2.10" 7 (A = -488MHz). For 
rubidium, k^ b = 8.05 x lOV^ 1 and T Rb = 6.1MHz. 

small cloud (koR <C 1), the scattering of light is isotropic, 
as expected from Rayleigh theory: only the first Mie 
mode (n = 0) is populated. However, for large many- 
particle spheres, Mie scattering turns out to be fun- 
damentally anisotropic, as many modes arc populated. 
Figflja) shows a typical momentum distribution of the 
excited atoms ^(p)! 2 and the associated scattering pat- 
tern |s c (k)| 2 of the light. Although most atoms recoil 
to k = kgz, a significant amount of light is scattered 
backward (k = — k^z) and acts as a seed for the MWSR 
instability. 



Note that the momentum distribution D e (p) = 
\ip e (p)\ 2 is concentrated along a circle with radius p = 
m c hko rather than ftfco (i.e., that 7n(p) reaches a max- 
imum for p = m c Hko). This is a signature of the 
Minkowski momentum for atomic recoil, that character- 
izes the momentum exchanged between light and matter 
in dielectric media 22 1. The blurring of the momentum 
wavefunction along the circle originates in the finite size 
of the cloud that creates a natural momentum spread 
dp ~ h/R. The ripples of the distribution are due to the 
sharp boundary of the cloud's density, yielding a Fourier 
transform with many secondary peaks. 

In optical dilute clouds, almost all of the light is 
scattered forward, leaving the incoming light almost 
untouched. Such a scattering pattern is displayed in 
Figfljb), where the intricate pattern of Mie scattering by 
an optically dense cloud has disappeared. In this limit, 
the excitation field is often approximated by the timed 
Dicke state /3(r) ~ /?TDse 4k °' r [HI, yet such an ansatz 
fails to predict any three-dimensional recoil pattern. 

After an atom absorbs a photon from the laser with a 
momentum kick m c ftko, the photon is reemited according 
to the Mie pattern s c (p/ft) oc ip e (p). Thus the atom will 
gain an extra momentum m c ftfc with a direction oppo- 
site to the emitted photon, and the momentum pattern 
of the ground state atoms after the scattering process is 
given by |^ e (m c ftko — p)| 2 - Experimentally, the column- 
integrated momentum distribution is observed in time-of- 
flight images. This leads to defining the projected distri- 
bution Dv(p x ,p z ) = J \tpg(p)\ 2 dp y . Such an integrated 
distribution is presented in FiglJJa). The atoms are ob- 
served to inhomogeneously fill a circle of radius m c hko. In 
particular, the part of the sphere around p = (3/2)ftk is 
weakly populated, which reflects the anisotropic nature 
of Mie scattering. 

Experimentally, the atomic recoil patterns are investi- 
gated by using the set-up of [l(|. After the interaction 
with the light, the 87 Rb atoms ballistically expand for 
tp = 20ms. Since the initial size of the cloud ~ 20/im 
is much smaller than its expanded size at the time of 
imaging, the density distribution of the expanded cloud 
is recorded by standard absorption imaging yielding the 
initial momentum distribution before the expansion ac- 
cording to p = mRbv/tF- 

The integrated momentum distribution observed ex- 
perimentally reproduces the features predicted by Mie 
scattering (see Figl^b)). The sphere of radius hko is 
filled with atoms, yet it exhibits a region where the proba- 
bility of the atomic recoil is very low around p z — 1.5ft/co. 
Moreover, it can be observed that a large number of 
atoms recoil around p = 2ftko, which is the signature 
of the backward emitted wave. The presence of preferred 
directions of emission is associated with the presence of 
Mie resonances [1(1]: the increase of emission in these 
directions can be interpreted as a sclf-Purccll enhance- 
ment [24j because the cloud acts as a cavity on itself to 
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modify its emission. This is in contrast to single-atom 
scattering, where light is scattered in random directions 
and not specifically backward and forward. 
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FIG. 2. (a) Integrated momentum distribution of the ground 
state, calculated from Eq |14l for a spherical homogeneous 
cloud of radius koR = 29.6, Ao = — 15GHz and with a refrac- 
tive index m c = 1.067 (a low-pass Gaussian filter was applied 
to attenuate the ripples due to the sharp boundaries of the 
homogeneous clouds, since they are irrelevant for comparison 
with the experiment) (b) Experimental integrated momen- 
tum distribution of the ground state for an ellipsoidal cloud 
of length koa z ~ 29.6, transverse radius ko<j± ~ 4.7, with 
N ~ 147000 atoms and a laser detuning Ao = — 15GHz and a 
20/xs laser pulse of 17mW. Its refractive index is m c « 1.067. 

It is important to note that Mie scattering is a three- 
dimensional process, while MWSR only develops along 
its most unstable direction. Mie scattering is a seeding 
process that emits light in many directions; including 
backwards, which is known to be the most unstable di- 
rection for MWSR in a cigar-shaped cloud illuminated 
along its main axis pjj . 

These two processes are illustrated in FigJ3J Initially, 
Mie scattering populates a sphere of radius p ~ m c hko- 
Then a MWSR instability develops, i.e., the grating in- 
duces light emission that, in turn, amplifies this atomic 
grating. The atoms observed in p rj — 2m c Sko and 
p w 4m c hk are not predicted by Mie scattering and 
can be explained only by the self-consistent matter-wave 
dynamics. 
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FIG. 3. Momentum distribution |V> 9 | 2 above the threshold of 
the MWSR. Experiment realized with an ellipsoidal cloud of 
transverse radius kocr± ~ 3.5, length koo z ~ 22, with N ~ 
156000 atoms and a laser detuning Ao = — 15GHz. 

In this paper we showed that Mie scattering induces 
a grating in the atomic distribution even below the 



threshold for the MWSR instability. Its signature is an 
anisotropic three-dimensional halo in the atomic momen- 
tum distribution. The atoms observed at p « 2Sko gen- 
erate a seeding wave for MWSR. Indeed, the matter wave 
modes at and 2fiko together form a density grating, at 
which subsequent light injected from the pump laser is 
Bragg-scatterered in a self-amplifying process. 

Note that quantum fluctuations and disorder may also 
give rise to backward emission [26| , and thus act as a seed 
for the matter- wave superradiance. These effects will be 
included in an extended work. 

Finally, it is interesting to remark that the off-axis 
emission of photons should be associated with higher 
modes (n 3> 1) that cor resp ond to photons with long 
lifetime within the cloud [27j- Thus, a time-resolved ob- 
servation of the off-axis atomic recoils should bear the 
signature of subradiance. 
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